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Abstract 

The transmission coefficient for vibrational waves crossing an abrupt junction between two thin elas- 
tic plates of different widths is calculated. These calculations are relevant to ballistic phonon thermal 
transport at low temperatures in mesoscopic systems and the Q for vibrations in mesoscopic oscilla- 
tors. Complete results are calculated in a simple scalar model of the elastic waves, and results for long 
wavelength modes are calculated using the full elasticity theory calculation. We suggest that thin plate 
elasticty theory provide a useful and tractable approximation to the full three dimensional geometry. 



1 Introduction 

The electronic properties of mesoscopic systems have been studied experimentally and theoretically for many 
years. More recently the behavior of other excitations, for example lattice vibrations (phonons) and spin 
degrees of freedom, have come under study in these systems. In this paper we present results relevant to the 
issues of heat transport by phonons and the dissipation of vibrational modes in mesoscopic systems. 

The interest in heat transport by phonons in mesoscopic systems arises because for easily fabricated devices 
the vifavelength of a typical thermal phonon becomes comparable to the dimension of the thermal pathway 
at accessible temperatures of order IK. Thus quantized thermal transport due to the discrete mode struc- 
ture of the thermal pathway should become evident. We and others H, || showed that this leads to a 
natural quantum unit of thermal conductance k^T/h similar to the role of e'^/h as a quantum of electrical 
conductance in one dimensional wires Q. This quantum unit of thermal conductance is predicted Q to 
be clearly observable at low enough temperatures where only the acoustic (w — > 0) vibrational modes are 
excited in the thermal pathway — the waveguide-like modes with nonzero frequency cutoffs at long wave- 
lengths are populated with exponentially small numbers: a universal thermal conductance is predicted equal 
to NATT^k^T/Sh with Na the number of acoustic modes (4 for a freely suspended beam of material, corre- 
sponding to a longitudinal mode, two bending modes and a torsional mode). These predictions were recently 
verified in beautiful experiments by Schwab et al. 

Vibrational modes in mesoscopic systems arc found to have anomalously low Q values, compared to larger 
systems of the same material [|[ 0, ^ ||, 0. At first sight, the dissipation might be expected to 
become smaller as the oscillator becomes smaller, since defects such as dislocations are eliminated when the 
size gets less than a typical defect separation. Thus the observation of lower values of Q was a surprise. 
In addition, unexpected dependencies on temperature pT| and magnetic field j^] have been observed and 
remain unexplained. 
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Both the possibihty of observing the universal thermal conductance and explanations for the Q of small 
resonators involve the properties of phonon excitations with a wavelength comparable to the system size. 
Here we investigate a particular issue relevant to both these questions, namely the coupling of vibrational 
modes across an abrupt junction between two blocks of the same material but with different dimensions. 



Cavity 



Bridge 



Reservoir 



Figure 1: Schematic of possible experimental geometry for the study of thermal transport and oscillations 
in mesoscopic systems. 



The geometry typical of a number of experiments is shown in Fig. The "bridge" is made of silicon, silicon 
nitride or gallium arsenide, is freely floating, and is of rectangular cross section. The bridge is connected 
to two larger blocks of the same semiconductor. In thermal conduction experiments the block at one end, 
called the cavity, is also freely floating (and is physically supported by four bridges) and is of the same 
thickness. The block at the other end provides both the mechanical support and a thermal reservoir. In 
recent experiments the dimensions were: thickness t = 200nm; width w = SOOnm; length L = 5/zm. In 
vibration experiments the bridge may be supported just at one end (a cantilever) or at both ends (a 
beam) 

An important question in both thermal transport and oscillator damping experiments is the coupling of the 
vibrational modes of the bridge to modes in the supports — how well the energy in a mode in the bridge is 
transmitted to the supports, and vice versa. We have previously introduced a simple scalar model of the 
elastic waves to study this question in the context of thermal transport In this paper we give a more 
realistic description of the vibrational modes. 

First, in section ^ we introduce an improved scalar model, using a better choice for the boundary conditions 
on the scalar field that provides a more realistic approximation to the waves in an elastic medium. The scalar 
model in its revised form provides a useful first guide to the expected behavior of the experimental system, and 
a simpler environment in which to develop intuition and methods of theoretical attack. With the scalar model 
we perform a complete calculation of the scattering of the waves at the abrupt junction between the bridge 
and the supports for all the modes and at all wave vectors. We use the resulting transmission coefficients 
to evaluate the effect of the abrupt junction on the thermal conductance, particularly the universal low 
temperature expression. In addition we introduce a simple method to calculate the transmission coefficients 
for the long wave length acoustic modes, and compare the results with the general results. In the full elastic 
calculation we will not be able to calculate the transmission coefficient for full range of modes, but will be 
restricted to this type of long wavelength calculation. 

Secondly we propose that the elasticity theory for a thin plate geometry provides a useful semiquantitative 
description of the experimental geometry. This is described in section ^. The full elasticity theory for the 
modes in the two dimensional, thin plate geometry is sufficiently tractable that a complete mode spectrum is 
readily calculated. On the other hand, a fully three dimensional elasticity theory can only be attacked purely 
numerically. The two dimensional theory reproduces many important features of a fully three dimensional 
elasticity theory, for example the mixing of bulk longitudinal and transverse waves by refiection at boundaries, 
the correct behavior of the dispersion relation at long wavelength and low frequency, including the "bending" 
modes with the unusual quadratic dispersion a; cx fc^ at long wavelengths, and regions of negative dispersion 
in the mode spectra. Thus the results should be more informative than the naive scalar model. The results 
should be accurate at sufficiently low temperatures where the modes with structure across the thickness are 
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frozen out. We use the thin plate model to investigate the mode structure in the beam, and the coupling of 
these modes to the supports, also treated as thin plates of the same thickness. 

Finally, in section^ we apply the results to the issues of heat transport and oscillations in mesoscopic systems. 
1.1 Heat Transport 



A thermal transport experiment is shown in Fig. (|^). Two thermal masses are connected by four thermal 
pathways of mesoscopic dimensions in which heat transport by phonons is the dominant mechanism. One 
of the thermal masses, which we call the cavity, is a freely suspended thin block of semiconductor, with 
resistive wires on the upper surface to act as heat source and thermometer. The four bridges act as the 
thermal pathway to the outside world, as well as mechanical supports. Conceptually, heat is added to the 
cavity by resistive heating, and the resulting temperature difference from the reservoir is measured, yielding, 
for small heating, the thermal conductance of the bridges. In practice issues such as the thermal contact 
between the electrons in the resistive heater and thermometer and the phonons, and other thermal pathways 
to the reservoir such as through the electrical contacts to the resistive heaters, have to be considered. In 
this paper we will focus on the ideal situation where the phonon thermal pathway of the bridge dominates 
the conductance. 

In mesoscopic systems it is easy to cool to temperatures where the transverse dimensions of the beam are 
comparable to or smaller than the typical phonon wavelength hc/ksT where ks is Boltzmann's constant, 
h is Planck's constant, c is a typical speed of sound in the material and T is the temperature. When this 
condition is satisfied the discreet mode structure of the thermal pathway becomes evident — the first level 
of the quantization of thermal transport. Angelescu et al. ||l| showed that the thermal conductance takes 
on a quasi-universal form, largely independent of the material and mode structure of the beam, on the 
assumption that the contact between the modes in the bridge and the cavity and reservoir can be considered 
ideal. An ideal contact implies that the right going phonon modes in the bridge in Fig. (|l|) are populated 
with a thermal distribution at the temperature of the cavity, and the left going modes at the temperature 
of the reservoir. 

In a thermal conductance measurement the cavity and reservoir are maintained at temperatures T + ST 
and T with temperature difference ST small compared to their mean temperature. If we first look at the 
transport by the right moving phonons, the energy flux is 



where k is the wave vector along the bridge, LUm{k) is the dispersion relation of the mth discrete mode of the 
bridge, and Vgm = dujra{k) / dk is the group velocity. Transforming the integral to an integral over frequencies 




Figure 2: Experimental geometry of Tighe et al. fl^ . 




(1) 
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yields an expression for the heat transport by right moving phonons 

1 

= — ^/ dwfu^r,-Xk)n{u^{k)), (2) 

where cj™ is the cutoff frequency of the mth mode, i.e. the lowest frequency at which this mode propagates. 
(We have assumed the mth mode propagates to arbitrarily large frequencies. If a particular mode only 
propagates over a finite band of frequencies, the upper limit of the integral will be replaced by w™'^''.) 
The key simplification in this result is that the group velocity factor is cancelled by the density of states in 
transforming from an integration over wave numbers to an integration over frequencies. 

For ideal coupling to the reservoirs the distribution function n{ujm{k)) for the right moving phonons in 
Eq. is evaluated as the Bose distribution at the cavity temperature T + 6T. The thermal conductance is 
given by subtracting the analogous expression H^~^ for the left moving phonons given by Eq. (^), but now 
with the distribution n {ujm (k)) given by the Bose distribution at the reservoir temperature T 

^^,^^H^-HT^ST)-Hi-HT)^ (3) 
ST^o ST 

Finally, introducing the scaled frequency variable x = huj/ksT gives the expression Q 

if= ^^/(fi^rnAsT) (4) 



where / is given by the integral 



Equation demonstrates the important result that the properties of the bridge only enter through the 
ratio of the mode cutoff frequencies to the temperature ficj^/fc^T. The quantity kgT/h plays the role of 
the quantum unit of thermal conductance, analogous to the quantum of electrical conductance e^/ft. for one 
dimensional wires. At very low temperatures the contribution to the thermal conductance by the modes 
with nonzero cutoff frequency will be exponentially small leaving a universal thermal conductance |^ 

K^N.^ (6) 

where Na is the number of "acoustic" modes, i.e. modes with frequency tending to zero at long wavelengths. 
Usually this will be four for the beam (two transverse bending modes, one longitudinal compressional mode, 
and a torsional mode). Note that there is no dependence on the bridge properties in this expression. 

More generally we cannot assume perfect coupling between the modes in the bridge and the cavity and 
reservoir. This can be taken into account, following the Landauer approach to electrical conductance Q], 
through a transmission coefficient for energy to be transported across the interfaces. For example for im- 
perfect contact at the cavity bridge interface we would find a thermal conductance (returning to unsealed 
quantities in the integral for clarity) 

where %n{uj) is the energy transmission coefficient from the mode m of the bridge at frequency uj into the 
cavity modes. Imperfect coupling at the bridge- reservoir junction, and elastic scattering due to imperfections 
in the bridge can be similarly included through a total transmission matrix as in the electron case. 

A central issue in predicting the thermal conductance is then to calculate the transmission coefficient 7^(a;). 
This is particularly important in the question of the observability of the universal conductance at low 
temperatures, since the scattering of the long wavelength phonons contributing to this quantity becomes 
strong — indeed for the abrupt junction in Fig. (|^), Tm{uj) ^ as w — > as we will see below. Although 
it is feasible in experiment to "smooth off" the corners, as indeed was done in the experiment of Schwab 
et al. j^, consideration of the worst case abrupt junction provides insight into the importance of geometric 
scattering. 
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1.2 Oscillator Q 



The Q of an oscillator is given by 



where E is the rate of energy loss from the mode at frequency lo containing energy E. If we consider the 
oscillations of a beam supported by two supports, or a cantilever with one support, and estimate the energy 
loss as the energy transmitted into the supports, we find for the mode n 

where Vg = duj / dk is the group velocity of a wave propagating in the beam, and L is the length of the beam. 
The exact evaluation of this quantity depends on the nature of the mode (longitudinal, bending etc.). For 
the longitudinal and torsional modes, which have a linear spectrum lo — ck, the frequency of the fundamental 
mode in a beam of length L is of order cn/L, the group velocity is Vg = c, and so 

Qn' - -■ (10) 

rnr 

For the bending waves with a quadratic spectrum a; cx fc^, the result is more complicated, but the Q values 
are similar, and tend to this form for large n, so we will use this expression as a fairly accurate general 
estimate. 

We see from Eq. ( p^ that good isolation of the mode 7^ — > is a criterion for high Q. In practice this 
expression for the dissipation may be an overestimate, since we are assuming that all the energy of the mode 
that enters the support either dissipates away, or propagates away to large distances so that the energy is not 
returned to the oscillations of the beam. If this is not the case, the transmission of energy into the support 
is only one part of the problem — we would also have to consider the behavior of the vibrational energy in 
the supports as well. 



2 Scalar Model 

As a simple model of the elasticity problem consider a single scalar field </>. This might represent, for example, 
the (scalar) "displacement" , and the (vector) "stress" E would then be proportional to V(/). We will suppose 
a two dimensional domain corresponding to the thin plate. This leads to a wave equation 

(u, 

with V~ the two dimensional Laplacian and c giving the speed of propagation of the wave. Stress free 
boundary conditions at the edges are then 

n-\7(f> = 0, (12) 

with n the normal to the edge. Note that this Neumann boundary condition allows the propagation of an 
acoustic mode {uj{k ^ 0) ^ 0) in the bridge as we expect for elastic waves, whereas Dirichlet boundary 
conditions do not. An example of an elastic system described by such a scalar model is a stretched membrane: 
(f) would then be the displacement normal to the membrane and n • S is the vertical force on a unit line in 
membrane normal to the direction n. 

The scalar problem is sufficiently simple that we can calculate the transmission across an abrupt junction 
such as the cavity to bridge junction in full detail. This allows us to gain insight into the more complicated 
elastic wave problem, and also allows us to illustrate and test approximation schemes that will be useful 
there. 

The model Eqs. (O, |l2|) was studied by Angelescu et al. Q using the mode matching method developed by 
Szafer and Stone |l4| for the analogous electron wave calculation. However Angelescu et al. implicitly used a 
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rather unnatural boundary condition (0 = 0) for the end of the cavity at the junction plane (although Eq. ( p^ ) 
was assumed everywhere else, i.e. on the edges of the beam and cavity parallel to the propagation direction). 
We briefly review this work, and explain how this boundary condition was introduced, and then treat the 
more natural case (Eq. ( p^ ) everywhere) using the same methods. This new treatment actually removes a 
weak logarithmic divergence found in the original treatment, and produces results at low frequencies that 
are more consistent with the results of the full elasticity treatment. 



2.1 Model of Angelescu et al. 



w 




Figure 3: Geometry for the calculation of the transmission coefficient Stress free boundary conditions are 
assumed on the edges as shown. Angelescu et al. used cf> = boundary conditions on the end AB and CD 
of the cavity. A better choice is to use stress free conditions here as well. 



Assume a simple two dimensional geometry consisting of a rectangular bridge of transverse dimensions w 
connected to a rectangular cavity of transverse dimension W, Fig. (^. In the general three dimensional 
case, if the cavity and bridge have the same thickness, there is no mixing of the z modes, and the problem 
separates into a set of two dimensional problems, one for each z mode. Here we will only consider the lowest 
mode with no structure across the z-direction, which is the only mode excited at low enough temperatures. 
Let Xaiy) ^^'^ Xmiy) be orthonormal transverse modes satisfying the stress free boundary conditions on the 
edges in the cavity and the bridge, respectively. (For clarity we will denote cavity mode indices by Greek 
letters, and bridge mode indices by Roman letters.) 

The solutions to the wave equation take the form 

0™(^,y,t) = X™(2/)e'('"-"*' (13) 

for the bridge, where the frequency of the mode lo is given by = + c^k^ with Um = rmrc/w the 
cutoff frequency of the mth bridge mode. The form is similar for the cavity modes with the cavity width W 
replacing the bridge width w. We will denote the frequency separation between bridge modes by A: 

A = LOm+l - UJm = TTc/w. (14) 

Consider a phonon incident on the interface from the cavity side (x < 0), in the mode a of the cavity, and 

(c) 

with longitudinal wave vector ka ■ The solutions in the cavity and bridge, including the reflected waves in 
the cavity and the transmitted waves in the bridge, are: 

= + r^pxfe~'''f^ cavity, (^5) 

(t> = Em iamXme*'''"^ bridge. 
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with r^p and tam reflection and transmission amplitudes to be determined. In the above equations, fc„ 

and k^p'' are the wa 
matching condition 



(c) 

and are the wave vectors of the transmitted and reflected waves respectively, given by the frequency 



u.^ = + ^ c'kl + ^ <^k^f + u.f\ (16) 

Note that the sums over m and (3 in Eq. ( p^ ) include evanescent waves (imaginary k or fc^'^^) although only 
the propagating modes will contribute to the energy transport. The field <f> and the longitudinal derivative 
d(j)l dx have to be matched in the medium at a; = 0, which leads to the equations: 



(c) I ('^) _ f 

Xa "T — 2^rn^O'^Xmi (^'7\ 

jL(c) (c)__y^ _,jL(c)^(c) _ f h. ^ ^ ' 



Equations for the reflection and transmission coefficient are extracted by integrating Eqs. (O) multiplied by 

(c) 

a transverse function, Xm or Xa , and using the orthogonality of the functions over the appropriate domain 
to extract relationships for the mode coefficients. We first multiply one of the equations with a cavity mode 

and integrate over the cavity width, making use of the orthonormality relation J dyx^a^X'p'^ — ^aP- 
this section we follow Angelescu et al. jl] and perform this operation on the first equation (i.e. the matching 
equation for the field (f>). It is at this stage that the boundary condition on the cavity field at the face x = 
for \y\ > w/2 is introduced. The replacement in the integration on the right hand side 



/W/2 i-w/2 
dyxi'U(x = 0)^ / dyxi'^Y.*amXm (18) 
-w/2 Jw/2 ™ 



rW/2 i-w/2 

dyxi'U(x = 0)^ / 

-w/2 Jw/2 

implicitly forces the boundary condition (j){x = 0) = for \y\ > i(;/2 



Multiplying the first equation in ( p^ ) by a cavity mode x^'^\ integrating over the cavity width, and using 
the orthonormality of the cavity modes, leads to 

''a/S = — <5a/3 + ^ ^QmOm/S, (19) 
m 

where amp is the overlap of cavity and bridge transverse functions 

amp = / dyxpXm- (20) 

J -w/2 

Equation ( p^ may now be plugged into the second equation in Eq. (p^, and the result is: 

'Xa ^ ~ ^ ] ^ ] tamO-mpk^p "'xj? ~ ^ ] tamkmXm, (21) 
m p m 

which, when integrated with Xm over the bridge width, yields 

n 

This is a system of equations that determine tam- bi Eq. ( p^ ) the kernel Amn is given by 

A-mn — ^ O-mpO-npf^p"^ ■ (23) 

These equations may be solved for the t's, and then the flux transmission probability from the wave vector 

(c) 

ka state of cavity mode a to bridge mode m is given by 



(24) 

ka 
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Now summing over all the cavity modes that are propagating at frequency uj leads to the "transport trans- 
mission coefficient" from the cavity to the mth mode Tm (uj) (for co > ajm) by 

{^) = ^ ^ "^am — ^ ^ |^Qm| • (25) 



(c) (c) "'O 



Wit 



ith km (cj) and {uj) given by Eq. (|l|). This also gives the energy transmission coefficient from the mth 



bridge mode to the cavity, by the usual reciprocity arguments. 

Equations ( |2^ , p5| ) involve sums over cavity modes. We may either evaluate the sums directly for a chosen 
value of the width ratio W/w, or may take the limit of a large cavity width W ^ oo when the sums are 
replaced by integrals. We calculate the matrix Amn for m,n < iV,„ax with N^^x some upper cutoff for the 
number of bridge modes retained and invert the A^max x -^max matrix system numerically to find tam and 
hence 7^ (lo). 

There is a simple approximation that provides an analytic form for the solution to Eq. ( p2| ) that is in 
reasonably good agreement with the exact solution. The approximation derives from three properties of ama ■ 
Firstly, ama = unless m and a have the same parity. In other words, even modes couple to even modes 
and odd modes to odd modes only. Secondly as a function of a, Gma is sharply peaked around a = mW/w, 
the width of the peak being of order W/w. And thirdly, must satisfy the completeness relation 

a 

The first two properties permit the key approximation, namely that Amn dmn (since the product of two 
functions peaked at different channels m, n is very small and Amn is rigorously zero when m and n are 
different parity modes). Then we only need the diagonal part of A 

Amm — ^ flm/jfc^ \ (27) 

which is simply a weighted average of the complex wave vector over the narrow range of reflected cavity 
modes for which am/3 is significant. (Note '^m/s = 1 by completeness). In this case Eq. ( |2^ ) separates 
into 

^ ^mOL — Amm^am ^~ ^ni^am (28) 

and then 



(29) 



The fiux transmission probability from cavity mode a to bridge mode m is given in this approximation by 

la,n-\tam\ ^(.) - (fc^ + ^^)2 + ' 

where Km = Re Amm and J,„ = Im Amm ■ The energy transmission coefficient is 



{km + Km)^ + Jm {km + Km)^ + Jm 



We now use the explicit form of the transverse modes to evaluate Tm (w). With the boundary condition 
d(t>/dy = at the y boundaries we have: 



/ 2 / rmry \ 

' — cos — - m even 
xUy) (32) 
'IsH'^) modd 
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a even 



Asin(H^) a odd 



(33) 



with m and a integers, with the special case Xoiu) 
calculated: 

ama = m, a not both even or odd 



l/y/w and Xo(?^) ^ l/VW- The ama are easily 




+ 



21V 2 



sin ( 



2W 2 ) 



2W 2 



m, a even 7^ 
m, a odd 7^ 



a even 7^ 



For large m, a (both even or both odd) we can approximate 



w sm 



(^ 

V 2 



2ty 



2 J 



2W 



2 



(34) 
(35) 

(36) 

(37) 

(38) 
(39) 

(40) 



The a„iQ are indeed sharply peaked as a function of cavity mode number a. The large m approximation is 
essentially identical to the result in the electronic case |TJ]. However the small, second term in the braces in 
Eq. ( p^ ) ignored in this approximation appears with the opposite sign in our application. This turns out to 
render the sum over the cavity modes (3 appearing in Eq. (^3|) weakly (logarithmically) divergent for large 
p. (Note that fc^'^^ is imaginary here, so this divergent contribution is to Im Amn and to the component Jm 
of the diagonal terms.) We must impose some upper cutoff to the sum to achieve finite results. Physically 
we might suppose such a cutoff may come from the breakdown of the sharp corner approximation at short 
enough scales, or ultimately, in a perfectly fabricated mesoscopic system, from the atomic nature of the 
material leading to a finite number of modes. 

Results for the transmission coefficient of the lowest bridge mode Ti^[io) are shown in Fig. (|^). There is 
only a weak dependence on the cavity sum cutoff. The dependence at small frequency fits well the expected 
[ p^ cubic dependence 7o(w) cx lu'^ . The results shown were calculated for the case of an infinitely wide 
cavity (sums over cavity modes replaced by integrals). Results for finite widths (e.g. W/w = 20) are very 
similar. The first 11 bridge modes (6 even modes) were retained in the matrix inversion for the results shown: 
increasing this number did not change the results significantly showing that Amn indeed decreases rapidly 
for increasing \m — n\. Note that 7o(a;) rapidly approaches unity as the frequency grows. There is a small 
(< 10%) decrease at the frequency uj — 2A where the second even bridge mode becomes propagating, and 
similar features of reducing size occur at subsequent integral multiples of 2 A. There is no coupling between 
even and odd modes for the symmetric geometry used. 

The comparison between the results from the full matrix inversion and the diagonal approximation is shown 
in Fig. (|^). The diagonal approximation gives results good to about 10%. However this comparison depends 
on the cavity wave number cutoff, since the diagonal approximation depends rather more strongly on this 
parameter than for the full results shown in Fig. (Q). Note that the feature at a; = 2 A due to the interaction 
between different bridge modes is absent in the diagonal approximation. 



2.2 Stress Free Ends 

2.2.1 Mode matching calculation 

We now redo the scalar analysis, enforcing the boundary condition h ■ Vcf) = on the end of the cavity a: = 0, 
\y\ > w/2. This corresponds to a stress free boundary everywhere. 
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Figure 4: Transmission coefficient coupling the lowest bridge mode to the cavity modes as a function of the 
reduced frequency of the mode w/A with A the splitting between bridge modes at zero wave vector. Curves 
for cavity wave number cutoff equal to 20, 40, 80, 100, 500 times ttw~^ show the weak dependence on this 
cutoff. 



The analysis proceeds as before up to Eq. ([Tsl). But now we first multiply the second equation (for the 
continuity of dfjj/dx) by x^p^ ^-^d integrate over the cavity width. This enforces the boundary condition on 
the cavity face 



(c) 

The orthogonality of the gives 



n.V(t> = for x = and |y| > w/2. (41) 



Use this equation to eliminate r^p from the first of Eq. ( p^ 

n 13 m 

and integrate with Xm over the bridge width to yield 

n 

where 

a 

Again we can solve the equation for ^ctm Eq. ( p4[ ) numerically or by using the diagonal approximation for 
Ajnn- It is easily seen that the extra inverse powers of ka ^ in render the sum over a convergent, unlike 
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Figure 5: Tranmssion coefficient for the lowest bridge mode in the Angelescu et al. scalar model calculated 
using the (exact) matrix diagonalization (points) and the diagonal approximation (crosses) for a cavity mode 
wave vector cutoff 407rit;~^. The diagonal approximation shows a somewhat stronger dependence on this 
cutoff than shown by the exact results, so that this comparison will vary as we change the cutoff assumption. 



the case for The energy transmission coefficient from the mth bridge mode remains given by Eq. ([25|). 
The diagonal approximation A„m — AmmSmn now leads to 

T (ui) ~ '±JVm"-m (46) 

with J„, = ReAmm and Km = lm^„„i. 

The result for Tq (ui) for the lowest bridge mode is shown in Fig. (|^) taking the limit of an infinite cavity 
width (the sums in Eq. ( ^5| ) etc. evaluated as integrals). Again the transmission coefficient grows rapidly, 
approaching close to unity as the frequency grows, for example reaching 0.9 by about uj ~ 0.5A. Note the 
important difference from the previous scalar model that the low frequency asymptotic behavior is linear, 
Tq (lo) (X UJ, rather than cubic as obtained there, and so Tq (w) approaches unity more rapidly than anticipated 
in that work. A small reduction in Tq (ui) (by about 3%) near u; ~ 2A and by smaller amounts at higher 
integral multiples of this frequency are apparent. An analysis of the curve in this region shows a square root 
dependence on a; — 2A, corresponding to a coupling in the full matrix calculation to the second bridge mode, 
and to the square root growth in Tz (lli) that occurs here. The diagonal approximation (not shown) gives 
results that are indistinguishable on the figure for < uj < 2A, but the small decrease above ui = 2A does 
not appear in this approximation. 

The inset in Fig. (|^) shows the comparison with results for a finite cavity width. The results are quite 
surprising, with resonance- like features occurring whenever the frequency passes through the cutoff frequency 
of a cavity mode. This can be traced to the inverse power of fc^'^' occurring in Eq. (pa) and Eq. (Eq). Since 



goes to zero as y uj — ujn'' these singularities are integrable, and the features disappear in the limit of 
infinite cavity width. Smoothing over the features (e.g. taking the average over bins between successive 
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(o/A 



Figure 6: Transmission coefficient coupling the lowest bridge mode to the cavity modes as a function of the 
reduced frequency w/A in the scalar model with stress free boundaries. The main graph is for an infinite 
cavity width. The inset shows the comparison with results for a finite cavity width {W — 20.217u'). 



uJn ) gives points that follow the smooth curve for the infinite cavity width closely. Integrating the effect of 
the transmission coefficient in the thermal conductance over the thermal factors will effectively perform this 
smoothing, so that the features will not be apparent in the thermal measurements. The sharp features are 
presumably also smoothed out if the junction is not perfectly abrupt. 

Results for Tn (to) for other values of n are shown in Fig. (^. 

It is now straightforward to calculate the thermal conductivity using Eq. (^. We focus on the low temper- 
ature limit where the universal result for K/T applies in the ideal limit. With no reduction in the thermal 
transport due to scattering at the junction the plot oi K/T as a function of temperature develops a plateau at 
low temperatures at the universal value T:'^kg/3h (see Fig. (^). In this regime the conductivity is dominated 
by the acoustic mode with — s- as fc ^ 0. Scattering at the junction reduces the transmission of the 
small uj,k waves, so that the value of K/T is reduced from the no-scattering value. As can be seen from 
Fig. ( p^ ) this reduction begins to occur as the temperature is lowered at about the same temperature at 
which the plateau in the ideal case begins to develop. This suggests that the plateau at the universal value of 
K/T will not be well developed for the abrupt junction. The full calculation using Eq. (|^) and the reduced 
transmission coefficients of all the modes (solid curve in Fig. (0)) shows that this is the case — including the 
effects of scattering at the abrupt junction yields a K/T curve that tends smoothly to zero as T ^ 0. This 
result clearly demonstrates the importance of using smooth junctions between the bridge and the reservoirs, 
such as was done in the experiments of Schwab et al. ||], if the universal value of K/T is to be apparent. 

2.2.2 Long- Wavelength Calculation 

Although we have performed the full calculation of the transmission coefficient for the scalar case, this will 
not be possible for the elasticity theory calculation. We therefore introduce a simple analytic technique for 
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2.5 



(o)-co^)/A 



Figure 7: Values of the transmission coefficients from the nth bridge mode to the cavity as a function of 
(iij — LJn)/^ with ujn the cutoff frequency for the nth mode and A the spacing between mode cutoff frequencies 
A = ujn+i ~ ^n- The long-dashed line shows the square root dependence near cutoff for the n ^ Q modes. 
The short-dashed line is the linear dependence T[){lo) ~ 27ra;/A expected for the lowest mode. 



establishing the low frequency limit of Tq (cj) that can be extended to the full elasticity description. 

The approach relies on the poor transmission To {lo) for a; — > (long wavelength waves are strongly 
affected by the abrupt junction) to treat the transmission perturbatively. First the bridge mode is calculated 
assuming perfect reflection i.e. isolated from the supports. A simple analysis shows that the appropriate 
boundary condition on the end of the bridge is zero displacement (rather than zero stress). If we now 
"reconnect" the bridge to the supports, the stress fields at the end face of the bridge act as radiation sources 
of waves into the cavity. The total power in these waves for unit incident amplitude in the bridge gives us 
the transmission coefficient. The total power radiated may be readily calculated by integrating the product 
of the stress source and resulting velocity across the end face of the bridge. 

To establish the appropriate boundary conditions in the To {lo) — > limit we first consider a finite cavity 
width i.e. an abrupt junction between a bridge of width w for x < and a cavity of finite width W ^ w 
for a; > in the limit w, fc — > 0. (It is simplest here to consider t he tra nsmission from bridge to cavity, and 



we have reversed the sense of the x coordinate compared to Sec. 2.2.1 .) For a; < we consider an incident 
wave of unit amplitude and a reflected wave of amplitude r 

cj) = e*'=^-^"* + re-'''''-"^\ x <0. (47) 
For a; > there is only the transmitted wave 

^ ^ ^^tkx-tu^t ^ x>0. (48) 

Here the wave numbers k are fixed by the dispersion relation uj = ck, which is the same on both sides of the 
junction for the acoustic mode. 

The reflection and transmission amplitudes can be calculated by a simple matching at a; = 0. (This is 
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Reduced Temperature 



Figure 8: Thermal conductivity divided by temperature reduced by the zero temperature universal value 
TT^kg/3h as a function of the reduced temperature kBT/HA: solid curve - full calculation including the 
transmission losses due to the abrupt junction for stress free face; long-dashed curve - contribution from 
the lowest (acoustic) mode, showing the reduction at low temperatures due to the scattering at the abrup 
junction; dash-dotted curve - ideal result from all modes with full transmission; dotted line - low temperature 
asymptotic slope predicted from the low frequency behavior of the transmission coefficient. The short-dashed 
curves shows the result including scattering of the calculation of Angelescu et al. [El. 



equivalent to a wave impedance calculation) . Matching the displacement field (f) gives 

1 + (49) 

and matching the total force gives 

w{l-r) = Wt. (50) 

Note that the force (the integrated stress) is conserved because there are no additional stresses on the cavity 
face for |?;| > w/2: this matching would not be appropriate for the boundary condition used by Angelescu 
et al. §. 

Thus we find 

W — w 2w , , 

r = , t = . 51 

W + w' w + W ^ ' 

These expressions are good for kW, kw <C 1. In this limit the matching conditions can be applied outside 
of the region close to the junction where the fields are perturbed from their asymptotic forms Eqs. ( |4^ , ^ ) 
but before the exponential phase factors of the wave propagation have significantly changed. For W ^ w 
these expressions reduce to r ~ —1, < ~ 0, i.e. strong mismatch and almost perfect reflection with a sign 
change. Note that at x = this implies ~ 0, i.e. zero displacement boundary conditions, together with 
d(t>/dx ~ 2ik for unit incident amplitude. We now use this result as the basis of the calculation of the 
transmission for k small but nonzero and the limit — > oo. 
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The transmission at /s 7^ is calculated as a radiation problem, namely via the power radiated by the end 
of the bridge into the cavity. The zeroth order approximation for the solution in the bridge is the perfect 
reflection result calculated above, giving the stress radiation source on the wall of the cavity 



s{y) 



d(j)/dx ~ 2ifce-*"* for \y\ < w/2 
for \y\ > w/2 



where the second line is just the stress free boundary condition 

d(j)/dx = for a; = 0, \y\ > w/2. 



(52) 



(53) 



The problem of the radiation due to a stress source on the boundary of an elastic half space is known as the 
Lamb problem in elasticity theory, and has been considered by many authors, for example reference [p^ . 
The radiation field can be calculated by standard Fourier transform techniques. The solution in the cavity 
for X > can be written 



1 



with 



_ r ^k^ - for Id < k 

" I Vfc' - for Id > k 



(54) 



(55) 



corresponding to propagation or exponential decay away from the interface. The transform 4>{C) is given by 
matching to the known d(j)/dx at x = yielding 



z#(C) 



s{y)e-<ydy. 



(56) 



The power radiated is the product of the stress d(t>/dx and the velocity d<j>/dt across the radiation source 



P = 



dt dx I 



where the <> denotes the time average. For the fields varying as e"*'^* this gives 



P = - Re 
2 



Inserting the Fourier expression for (j)(x = 0) 



w/2 



w/2 



4>{v)s*{y)dy 



x=0 



F = — Im 

47r 



dCm / dys*{y)e'^y 



47r 



Im 



dC- 



iq 



— 00 
00 



3{y)e-<ydy 



(57) 

(58) 

(59) 
(60) 



Since the transverse wave vector ^ is limited hy C, < k for the integrand to be imaginary (corresponding to 
q real, i.e. propagating waves in the cavity), and s{y) is nonzero only for \y\ < w/2 an expansion in kw <^ 1 
is given by the expansion 



s(y)e ^''^dy = sq - iCsi + • 



where Si are successive moments of the source 



So 



Si 



s{y)dy, 
s{y)ydy. 



(61) 

(62) 
(63) 
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We will need only the leading order term, i.e. sq if the source is parity symmetric in y, si if the source is 
antisymmetric. In particular for the present case sq — 2ikw so that 



Normalizing by the power in the incident wave ■^wujk gives the transmission coefficient 

4kw 1 



%{lo) = / ^==dC (65) 

The integral expression for Tq can be easily understood physically as integrating over the power radiated into 
the continuum of waves, traveling at the wave propagation speed c, and propagating at all angles into the 
half space. For the scalar model there is a single type of propagating wave. With the full elasticity theory 
we will see a similar result, but with a number of propagating waves contributing to the power radiated. 

Performing the integral gives 

To{uj) = 2kw = 2ujw/c for kw < 1. (66) 
This result confirms the linear dependence on uj for small u, as shown in Fig. (FtI). 



3 Thin Plate Theory 

A useful model of the mesoscopic geometry, that is more tractable than a fully three dimensional elasticity 
calculation, is to assume a thin plate geometry. Thus we take the elastic structure to be carved from a thin 
plate of uniform thickness d, which is taken to be small with respect to the other dimensions and also with 
respect to the wavelength of the elastic waves. In this model the mode frequency cutoffs at fc = 0, important 
in the thermal conductance, can be readily calculated (for most of the modes they are given by simple 
analytic expressions). In addition, although the mode structure is quite complicated, involving the mixing 
of longitudinal and transverse components by reflection at the edges, explicit expressions for the transverse 
structure of the modes in terms of a finite sum of sinusoidal functions can be written down. Compare this 
with a full three dimensional analysis, where a finite dimensional representation of the modes is not possible. 
In addition to the scattering at an abrupt junction pursued here, the thin plate limit will permit an analysis 
of phenomena such as the scattering of the waves of surface and bulk imperfections, an issue we hope to 
pursue in future work. 

We first review the general elastic theory for waves in a thin plate, and confirm the expected dispersion 
relations for the acoustic modes of a rectangular beam (i.e. the "bridge") in the long wavelength limit, and 
then use the equations to study the coupling of long wavelength modes across the abrupt bridge-support 
junctions. 



3.1 Review of Elastic Theory and Modes 

The elasticity of an isotropic solid is summarized by the relationship between the stress tensor T and the 
strain tensor 

T,, = -ife% -2MSy. (67) 
Here Q is the dilation and S is the shear strain 

dx dy dz 
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with u(x) the displacement and K and /i elastic constants. 



For a thin plate of thickness d in the xy plane, linear elasticity theory can be separated into equations for the 
normal displacement Uz = w{x,y) of the plate and for the in-plane displacements averaged over the depth 
u{x, y) = (it, v) with u = {ux{x, y, z)) ^ and v = {uy{x, y, z)) ^, all functions of just two spatial variables. This 
is done by assuming, for in plane wave vectors k such that fcc? ^ 1, that the stresses in the vertical direction 
Sj;j-, which must be zero at the nearby stress free top and bottom surfaces, may be put to zero everywhere. 
This allows the variation of the strains across the thickness of the plates to be eliminated in terms of the 
variables u, v, w. For example setting T^z to zero gives 

duz _ K~lfi /du dv\ 

dz K+lti\dx^ dyj- ^ ' 

Thus the modes separate into modes with in-plane polarizations, and modes with polarizations normal to the 
plane (flexural modes). The full development can be found in any standard text on elasticity, for example 
Landau and Lifshitz lITstl or Graff IITtI. Here we collect the main results. 



3.1.1 In-plane polarization 



Using relationships such as Eq. (69) leads to an effective two dimensional elasticity theory for displacements 
in the plane, summarized by the stress-strain relationship 

^ _^0(2)^^^. _ 2/2S]g\ (70) 

where the indices i,j now run only over x and y, 0'^^ and S'^^ are the two dimensional dilation and shear 
strain tensor 

9(2) = ^ + ?^ (71a) 
ox oy 

and the effective two dimensional elastic constants are 

p. = ^J, (72a) 

The elastic waves with in-plane polarization are then given by the equation of motion 

= -V. • T(^) (73) 
= A'Vi (V_L • u) + /iViu. (74) 

with Vj_ the horizontal gradient operator. In a horizontally infinite sheet there are longitudinal and transverse 
waves with speeds 

p P (A + 1^) 

4 = ^ = ^. (75b) 
P P 

Alternatively, introducing the Young's modulus E and Poisson ration a (with — 1 < cr < 1/2) so that 
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we have 



Rather than using the Poisson ratio directly it is convenient to introduce the parameter for the ratio of wave 
speeds 



r='-^ = J-^>l. (78) 
ct V 1 — o" 

In terms of the component displacements u, v and the parameter r the elastic wave equation can be written 

1 d'^u 2^^" ^ + ( ^ 1) (79 
c|, dt^ dx^ dy"^ dxdy ' 

c\ dt^ dx^ dy^ dxdy' 

A result that we will find useful later is for the component displacements in the solutions for the propagating 
waves 

f uo{l,—kTx/kTy)e^^^^'^~"*^ transverse wave , . 

\ uo{l,kLy/kLx)e'^^^-''-'^*^ longitudinal wave ' ^ ' 

We note for completeness that in a three dimensional elastic medium, waves polarized in the xy plane having 
no z dependence are also described by an effective 2-dimensional elasticity theory as in Eq. @) but now 
with effective elastic constants K = K + ^fj,, p, = fi so that c|/c^ — 2(1 — o'))/(l — 2a). The difference in 
the effective elastic constants in the two cases arises because in the three dimensional medium restricted to 
no z-dependence, there can be no expansion in the z-direction to relieve the stresses set up by the strain 
in the xy plane. In many elasticity textbooks, the "two-dimensional elasticity" and "thin-plate theory" 
discussed correspond to this case. We can relate these results to the thin plate geometry we are considering 
by appropriately transforming the effective elastic constants. 

For a finite plate we must apply stress free boundary conditions at the side edges 

n-T(2^=0 (81) 

with h the normal to the edge. For waves propagating in the x direction in a long finite plate of width w 
the no-stress boundary condition at the edges of the plate at y = ±w/2 are therefore 

nf'O'l-ir''^ + 2)^)^ (82b) 



The boundary conditions have the effect of reflecting incident longitudinal waves into both longitudinal 
and transverse reflected waves, so that these waves become coupled in the finite geometry, leading to a 
complicated dispersion relationship. The solutions propagating in the x direction decouple into either even 
or odd signature with respect to y reflection, and take the form (using Eqs. (|80|)) 



cos(xTy) + a^l^ cos{xLy)\ e'C^^-*) (83a) 
= \-ik/xT)aP sin(xTy) + {iXL/k)a^l^ ^HxLy)] e^^'^-""*) (83b) 
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and 



{iklxT)a'T cosixTV) - {ixLlk)a^L cos(xLy) 
where Xt.l are given through the dispersion relation 



^i{kx—ijjt) 



(84a) 
(84b) 

(85) 



and we have defined K = lo/ct, the wave number of a transverse wave at the frequency uj in an infinite 
plate. The values of xt,l rnay be real or imaginary. Note that each wave combines both shear (xt) and 
compressional (xl) distortions, which are mixed by the reflection of the plane waves off the edges. 

The amplitudes a'"j)j^°'' must be adjusted to satisfy the boundary conditions aX y = ±w/2, Eqs. (|8^). 
Substituting into these conditions leads to a system of two homogeneous equations for ■* , a^^' (and two for 
CL'i/') and so a consistency condition that leads to a transcendental equation for w for each fc, known as 



the Rayleigh-Lamb equations]^. 



% 20 



o 




Figure 9: Dispersion relation of in-plane polarized modes for the thin plate. The value of the Poisson ratio 
is (7 = 0.33. 



For the even signature modes the transcendental equation is 

(fc2 _ tan ^ + ^exTXL tan ^ = 0. (86) 
For the odd signature modes the transcendental equation is 

A, 2 , XTW 2 2n2, XLW 

4k XTXLiari^- + {k - Xt) tan 0. (87) 

For a given value of the speed ratio r these equations can be solved numerically for uj{k). The spectrum for 
the value of r = -s/S corresponding to GaAs (cr ~ 1/3) is shown in Fig. ([l0|). 

'^The Rayleigh-Lamb equations usually occur in the somewhat different context of the analysis of the modes in a plate which 
is considered infinite in the xy plane, with propagation in the x direction and no dependence on the y coordinate. The wave 
numbers XT,L then give the variation across the thickness of the plate. The equations take the same form, with the wave 

speed ratio r given by the expression for two dimensional elasticity theory, r = 2{1 — a)) / {1 — 2a). For example Rego and 
Kirczenow plot dispersion curves for these modes across the thickness of an infinite plate. 
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At large values of k the slopes of the curves (except the lowest two) asymptote to cl or ct corresponding 
to the freely propagating waves in the plate. The lowest two modes asymptote to a slope cs less than both 
Ct and cl- In this case for large k we have tan{xT,Lw/2) ^ i in Eqs. ( |86| , p7| ) so that the slope cs = rscr is 
given by the solution of 

4^1^^1371/^ = (2 -r|)2. (88) 

This is an edge wave analogous to the Rayleigh wave on the surface of a three dimensional slab of material. 
For r = \/3 Eq. gives Cs = — 2/\/3ct- 

The values of the finite-frequency intercepts of the "waveguide" modes for ^ can also be calculated 
analytically. For the even mode intercepts Eq. ( ^6| ) is satisfied at fc ^ by tan(xTif/2) = or tan(xLU'/2) 
CO. Similarly Eq. ( p7| ) is satisfied by tan(xiw/2) = or tan(xTw/2) oo. Thus the zero wave number 
intercepts are given by the simple expressions for transverse and longitudinal wave propagation 

Lo';^^ — mrcT/w, u)^^ = rmiCT/w — n-KCi^jw. (89) 

(These simple results hold because for fc = there is no interconversion of longitudinal and transverse waves 
on reflection at the edges.) The shapes of the curves between A: = and the large k asymptotes are quite 
complicated, with various mode crossings and regions of anomalous dispersion diu/dk < for example. 




Figure 10: Dispersion curves for the long wavelength inplane modes of a thin plate beam of width w for 
Poisson ratio 1/3 (wave speed ratio VS)- Solid lines: numerical results; long dashed line: small k linear 
dispersion for compression mode; long-short-short dashed: linear dispersion for longitudinal wave in infinite 
plate; short dashed: linear dispersion for edge mode; long-short dashed: quadratic dispersion for small k 
bending mode. Note the anomalous dispersion of the fourth mode at small k. 

We are particularly interested in the long wavelength modes fc — > 0. The dispersion relation in this limit can 
be found by Taylor expansion of the tan functions in Eqs. ( ^6|j8^ ). The even mode tends to 

uj 2\J\ - r-^crk = ^/Ej~pk- (90) 

This agrees with the usual expression for the stretching mode of a rod. This is not the same as the dispersion 
for the bulk longitudinal modes in the thin plate, uj — CLk, but the speeds are quite close for cr = 0.33. On 
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the other hand the odd mode gives a quadratic dispersion, characteristic of bending modes of beams, 



y^^^^Tk^- (91) 

(This is given by expanding the tan functions up to cubic order.) Rod bending theory gives the expression 
uj/k'^ = yjEIjpA with / the areal moment of inertia about the midhne, and A the cross section area. For 
the rectangular beam we have I jA = /12, and using (1 — r~^)c|, — E/2p shows the correspondence. 



3.1.2 Flexural modes 



The flexural modes are most easily derived by using relationships such as Eq. (69) to derive an expression 
for the energy of transverse displacements Pq] 



where 



2(1 



D 




\dxdy J dx^ dy^ 



dxdy 



12(1 -(t2) 



(92) 



(93) 



is the flexural rigidity of the plate of thickness d. The equation of motion and boundary conditions are given 
by taking the variation of the energy with respect to displacements w(x,y). For a region with rectangular 
boundaries at x = ±a/2 and y = ±6/2 the variation is 



SF = D 



J J V^cj 6w dxdy 



dy I^Sw 



w 



dx^ 



dx I Sw 



dxdy"^ 



d6w 
dx 



dy^ 



(2- 



dx"^ 



dx^dy 



dy'^ 

dSw 
dy 



x=a/2 
x=-a/2 

a 



dy'^ 



dx^ 



y=b/2 



y=-b/2 



(94) 



The first term gives the effective force per unit area on the plate, and hence the equation of motion 



pd- 



9i2 



DWlw = 0. 



(95) 



The last two terms on the right hand side are boundary terms given by integration by parts. Physically they 
give the work done at the boundaries by the vertical force per length of boundary Vi against the vertical 
displacement dw and by the moment per unit length Mi against the angular displacement of the plate^ 
— ViW. Thus at X = a/2 we have 



and a.t y — 6/2 



Vx 



Mr. 



D 



dx \ dx"^ 
(d'^w 

K-dx^^' 



+ (2- 

d'^w 
dy"^ 



dy^ 



d'^w 



M, 



y 



-D 



dy \ dy"^ 
d'^w 
dy'^ 



+ {2 -a) 



d'^w 
dx'^ 



a 



d'^w 
dx"^ 



(96a) 
(96b) 

(97a) 
(97b) 



^ The sign convention for the moments Mx , My is that Mi is positive if it tends to produce compression in the negative z 
side of the plate. Thia^ angular displacements 8x,Qy are defined with the same convention. This is the usual definition in the 



elasticity literature |l8|. 
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For free edges, these quantities must be set to zero. In addition to the force per unit length there are also 
point forces localized at the corners e.g. ai x — a/2,y ~ ±6/2 fia] 



Fc ^ ±2D(1 - cr) 



dxdy' 



(98) 



These must be included when we are calculating the total force acting across the width of the beam, for 
example: 



F 



w/2 



V,dy + F,{w/2)+F,{-w/2) 

w/2 

= — / dy 

w/2 



dx 



(99) 
(100) 



and the latter equality shows the consistency with the macroscopic equation for the rotational equilibrium 
of the beam (see section 3.2.2| below). 

We now calculate the modes propagating in the x direction the bridge of width w. 

Again the modes have either even or odd signature with respect to y reflections. Since the wave equation 
is fourth order in the spatial derivatives, for each frequency lo there are two even or odd components. The 
solutions to the wave equation are (even) 



and (odd) 



where 



and we have written 



o'^ cosh(x+y) + a}- cosh(x-j/) 



a^°^ smh{x+y) + aL°^ sinh(x-2/) 



^i{kx—Lot) 



sJpdjDuj = K^. 



(101) 



(102) 



(103) 



(104) 



Again the dispersion cij(fc) and the ratio of amplitudes a+/a_ are determined by the requirement of consis- 
tency with the boundary conditions at the edges y — :kw/2 Eq. (|97|). This gives for the even modes 



(105) 



[K^ + (1 - cf)k^?x- tanh(x-u;/2) = [K'^ - (1 - cf)k^?X+ tanh(x+u)/2) 
and for the odd modes 

[K"^ + (1 - a)k'^fx- coth(x-u;/2) = [K"^ - (1 - cr)fcTx+ coth(x+u;/2)- (106) 



For LO ^ Q we can expand the hyperbolic functions and solve algebraic equations to determine the dispersion 
curve. 

For the even mode this gives K'^ /k^ — \p\~~a^ yielding the quadratic dispersion of the beam bending mode 

E 



pd 




dk' 



(107) 



agreeing with the expression from simple rod theory. We can follow this mode to large k where we find again 
a quadratic dispersion but with a different slope 



'l-3cr + 2^1-2cr(l-CT) 

(l-a)2(3 + fT) ■ 



(108) 
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Figure 11: Dispersion relation of the flexural modes of a thin plate with Poisson ratio a ~ 0.33. Note 
that the modes have approximately a quadratic wave number dependence (jj ~ o;^ + akF' with a ~ D / pd, 
corresponding to the "bulk" flexural wave. The two lower modes have a different form: one mode (the torsion 
mode) has a linear dispersion relation at low frequencies uj (x k and both modes asymptote to uj — (3k'^ at 
large k with (3 < a corresponding to an edge wave. It turns out that a c:i (3 for the particular value of a 
used here. 



This is again an edge wave (now an edge bending wave). The intersections of the higher modes with the 
frequency axis are given by k"^ ^ 0, x_ — *■ iK, x+ ^ K so that Eq. (105) reduces to 



ia.n{Kw/2) = tanh(if u;/2). 



(109) 



The solutions are well approximated by Kw = (3/2 + 2n)7r (i.e. u) = •\/D/pd(3/2+2n)^7r^/w^) for n = 0, 1 
(with an error of less than 1/2% for the worst case n = 0). 



For the odd mode Eq. (106) gives for fc — > the dispersion relation for the torsion mode 

t^^,, 2v/6(l-a)--2,/^-fc, 
pa w \ p w 



(110) 



which agrees with the usual result calculated in elastic rod theory. The large k asymptote of this mode is 
the same as Eq. (108). The intersections of the higher odd modes with the frequency axis are from Eq. (IOC) 
given by 



tan(ifw/2) = tanh(A'w/2). 



(Ill) 



The solutions are well approximated by Kw = (5/2 + 2n)7r (i.e. uj ~ yjD/ pd{^/2 + 2n)^7r^) for n = 0, 1 
Combining the even and odd modes, the zero wave number frequency intercepts can be written 



UJn — Ct- 



-.{-+nfn^d/w'^. 



The dispersion curves for ct = 1/3 corresponding to GaAs are shown in Figs. (pT|,|l2|). 



(112) 
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Scaled Wavenumber kw 



Figure 12: Dispersion relation as in Fig. ( |l l| ) but at small k. Note that the abscissa is proportional to k 
in this plot. Solid lines: numerical results; long-dashed line: expected quadratic dispersion for rod-bending 
mode; short-dashed line: expected linear dispersion for rod-torsion mode. 



3.2 Transmission coefficient in the infinite wavelength hmit 

The transmission coefficients for the acoustic modes in the long wavelength limit and for finite cavity width 
can be readily calculated by the wave matching methods as for the scalar waves. We investigate the trans- 
mission of very long wavelength modes at the abrupt junction between the bridge {x < 0) of width w and the 
cavity {x > 0) of finite width W. It is easiest to evaluate the wave fields by simple macroscopic arguments. 
As well as providing the long wavelength limit of the transmission coefficients, these results also provide the 



basis for calculating the leading order finite k corrections, following the same methods as in Sec. 2.2.2 



3.2.1 Compression modes 

The compressional (extension) modes for fc — > are given by the simple one dimensional calculation: 
where = E/ p. 

For a; < we have incident and reflected waves 

M= (e^'=^-Kre-*'=^)e-*"* (114) 

and for a; > a single transmitted wave 

u = ie*'=^e-*"* (115) 

where uj/k = ce- We match the displacement u and the total force either side of the interface 

l + r = t, (116) 
ikw{l - r) = ikWt. (117) 

Note that the force matching requires that the end surface of the cavity be stress free. This gives 

1 h-1 , , 

t= r, r^-- . (118) 

1+h h+l ^ ' 
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where h is the width ratio h = W/w. In the hmit ft, — )■ oo we find r — > — 1, i — > 0, i.e. perfect reflection with 
a sign change of the displacement. This implies that at the junction x ~ we have m ~ and the stress 
Edu/dx ~ 2ikEe-''^K 



3.2.2 Bending modes 



At long wavelengths the bending modes are given by equations for the total force F and the total moment 
M on each cross-section (see The moment from opposite forces on each end of a small element of the 

beam must cancel the net moment from the forces on the faces 



F : 



dM 
dx 



(119) 



(Since the moments scale as the length of the element 5x, whereas the moment of inertia scales as 5x^ , there 
is no inertial term in this equation). The net force on an element gives its acceleration 



pA 



dF 
dx ' 



(120) 



where u is the bending displacement and A is the cross section area. The moment M is given by the moment 
of the tensile stresses due to the extension and compression of the beam from its curvature 



M 



dx^ ' 



(121) 



where / is the areal moment of inertia of the beam around the midline normal to the displacement. These 
equations together give the equation of motion 



9*2 



Eld^u 
pA dx^ ~ ■ 



(122) 



The dispersion relation is quadratic, lj = ak^ with a — \J EI / pA. For a rectangular beam I /A = d? /12 
with d the thickness of the beam in the direction of the displacement. As we saw in the previous section, 
this reproduces the long wavelength limit of the acoustic bending modes calculated in thin plate theory. 



At a frequency cj, as well as the propagating modes at wave number i-y/cj/a there are also evanescent modes 
with decay rate zti/w/a: the modes localized at the junction and decaying to ±oo must be included in the 
mode transmission problem. 

For the bending mode with displacement normal to the plane the dispersion relation is the same in the bridge 
and cavity. Thus for an incident wave g^i'^^-'^t) in the bridge with k = ^JuTfa. we have 



te 



ikx 



re 
f he 



■ikx 
— kx 



ae 



kx 



X <o 
x>0 



(123) 



At the junction we require continuity of the displacement u, the rotation angle du/dx, the total moment 
—Eld^u/dx^ and the total force —Eld^u/dx^. This gives the matrix equation 



(124) 



where h — Id h — W/w is the ratio of the appropriate moment of inertia for the cavity and bridge. In the 
limit of large h the solution is easily found to be 



1 


-1 1 


-1 




r 




' -1 " 


1 


1 i 


i 




t 




1 


1 


-h -1 


h 




a 




-1 


1 


h —i 






b 




1 



r = I, 

t ^ i/h, 

a = -(l + i), 

6 = 2(1 - i)/h. 



(125a) 
(125b) 
(125c) 
(125d) 
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For the bending mode with displacements in the plane of the plate the dispersion relation is different in 
bridge and cavity. For frequency lu the wave numbers in the bridge and cavity are kb, kc with 



kjkb = K/a,)i/2 = (w/Wy/\ 
On the other hand the ratio of the moments of inertia is 

Ic/h = {W/wf. 



(126) 



(127) 



Thus the continuity of the displacement u, the rotation angle du/dx, the total moment —Eld^u/dx^ and 
the total force —Eld^u/dx^. This gives the matrix equation (writing h = ^JW/w) 



(128) 



1 


-1 


1 -1 


1 




i ih~^ 


1 




-1 


1 




-i -ih^ 





r 




' -1 " 




t 




1 




a 




-1 




b 




1 



In the limit of large h the solution is 

r = i, 
t ^ 2/P, 
a = -(1 + i), 
h = 2/P. 

With these expression in both cases we have at a; = to leading order in /i^^ or h" 

w = 0, 
du/dx — 0, 
d'^u/dx^ = -2(1 + i), 
d^u/dx^ = -2(1 + i), 



(129a) 
(129b) 
(129c) 
(129d) 



(130a) 
(130b) 
(130c) 
(130d) 



so that the displacements, u and the angle du/dx, tend to zero, but the corresponding stresses are large. 
Note that although the force and moment are out of phase for the single wave e^('^^~'^*), because of the 
evanescent waves near the junction they are in phase at the junction plane a: = 0. Equations (130) become 
the zeroth order input for the radiation calculation at non-zero wave number. 



3.2.3 Torsion modes 

The long wavelength limit of torsion waves is described by the one dimensional wave equation giving the 
angular acceleration in terms of the torque 



dP dx"^ 



(131) 



Here C is the torsional rigidity giving the torque on a section due to the twist r = dO/dx 

torque = Ct. (132) 

It is given by |l5j 

C = 4:^iJ xdxdy (133) 

where x satisfies the equation in the cross section 

Vix = -1 (134) 
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and X = on the boundaries. The form of the solution % is the same as the profile of the flow of a viscous 
fluid through the section and C is then proportional to the integrated flux. For the thin plate geometry 
with thickness d and width w ^ d, the value of C is ^fid^w and value of / is j^dw^, so that the ratio of 
propagation speeds in bridge and cavity is again the width ratio h — W/w. 

It is interesting to evaluate the stress distribution for the thin plate. The stresses in the (y, z) section are 
given in terms of x by 

ayx = 2^Tdx/dz (135a) 
(Jzx = -2fj,Tdx/dy. (135b) 
The solution for x is analogous to Poiseuille flow, so that 

X=^{d^-4z^) (136) 

except within a distance S ^ d from the side wall, where x must decrease to zero. Thus there is a distributed 
stress acting in the y direction 

~ -2^rz (137) 
and a stress in the z direction that is effectively localized (within a distance d) at the edge 

~ ^^^T{d^ - Az^)S{y - w/2). (138) 

This localized stress corresponds to the corner forces Eq. (^) in the thin plate theory. 

For an incident torsion wave 6 = ^^C^bx-uji) ^j-^g bridge for a; < we have incident and reflected waves 

and for a; > a single transmitted wave 

9 = te^'^'^^e"*"* (140) 

where kc/kf, — W/w. The matching of the angular displacement 9 and torque Cd9/dx for an incident wave 
gives 

l + r = t (141a) 
ikbw{l - r) = ikcWt (141b) 

so that 

with h — W/w. In the limit h ^ oo we find r —1, t 0, implying 9{x = 0) ~ and the stress 
Cd9/dx ~ 2ikCe-^'^K 

3.3 Transmission coefficient for small k 

In this section we calculate the small wave vector asymptotic limit of the transmission coefficient from the 



four acoustic modes of the beam into the cavity. The method follows that of section 2.2.2. Thus we calculate 
the radiation from oscillating stresses s(?/)e~"^* on the edge of the cavity. The stresses are calculated as 
the stresses arising on the ends of the bridge for zero displacement boundary conditions, as follows from the 
analyses in the previous section of the modes at infinite wavelengths coupling into a cavity of finite width. 
For the long wavelength value of the transmission coefhcients, only the radiation by the integrated stress 
J s{y)dy for the even parity modes, or the integrated first moment J s{y)ydy for the odd parity modes, 
is needed. The Lamb problem of the radiation from surface sources into an elastic half space has been 
much studied in the literature, for example see p^ ], and the reader is referred there for a more exhaustive 
discussion of this aspect of the calculation. The details of the calculations are quite complicated, and the 
reader may choose to skip these sections and refer to the discussion of the results tabulated in section ^ and 
the summary in table |l| there. 
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3.3.1 In-plane Compression 



For a compressional wave in tlie bridge of unit incident amplitude in the displacement, the oscillating end 
of the bridge acts as a stress source on the cavity face of amplitude 2iEk over the source region \y\ < w/2, 



embedded in the otherwise stress free line a; = 0. The solutions to the wave equations (79) in the cavity can 
be written (cf. Eq. (|^) 

u=— [aT(C)e'''^'' +aL(C)e'''^1e*^^e-*"*dC (143a) 
y^TT [-(fcT/C)«T(C)e^'"" + (C/A:L)aL(C)e^"-^e<^e-*"*dC (143b) 
where kx and are the z-components of the wave vectors of the transverse and longitudinal components 



ly/C^ - \C\>K ' I iy/C~Kyr^ \C\>K/r 



Id <K ^ ^ K^/r^ - Id < K/r ^^^^^ 



with — Lo^ /c^. (The signs chosen correspond to waves propagating away or exponentially decaying.) 
The amplitudes aL,T are fixed by matching to the normal and tangential stress sources for |y| < 1/2 

M^2?^ + (r2-2)|^) =!]„, (145a) 



dx dy 

In the present case I]„ = 2iEke~'^'^* for unit incident wave amplitude and — 0. Both components of the 
stress are zero for |y| > w/2. Taking the leading order expansion in Km for the Fourier transform of the 
source stress as in the scalar calculation, section [2.2.2| , gives for the Fourier components 

r'^{ikTaT + ik^aL) + (r^ — 2){—ikTaT + iQ^aL/kL) — Y^nw/ ^ (146a) 
iC{aT + aL) + {-ik^ar/C, + iCa^) = (146b) 
These equations are readily solved for a^^T from which we can calculate u{x — 0) for unit driving stress 

!fc^ = -!:!^i r^de (147) 

where ^ = rC^/K is the y-wave vector scaled by the longitudinal wave number K/r = lu/cl, 

FoiO = - r^f + ^e^T{^)^L{0 (148) 

and 



(This result is analogous to Eq. (123) of Miller and Pursey who calculate the average displacement 
at the aperture per unit oscillating stress for a line on a three dimensional half space. Indeed we can use 
their result if we express it in terms of the ratio of wave speeds, and the elastic constant /i which retains 
its significance unchanged between the two geometries. The translation from their (MP) notation is then 
Ux,MP {^/cl)ux, ciMP {(^/cl)w/2, iiMP r, C4^a.mp — > 1^- Note carefully that the usage of "/i" 
is different in their work and ours. We have also taken the leading order term in ujw/cl by making the 
replacement e*^^ 1 for |y| < w/2.) 

For unit incident wave in the beam the longitudinal stress at the aperture is E„ — 2iEke~^'^* and the power 
radiated is ^w'Re(— uS*). In the incident wave the stress is Sji, = iEke~^^*' and the velocity ii^ is 
so that the incident power is ^wEujk. The transmission coefficient is therefore 



-toje 



%{u)^^{kw) 



(150) 
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where we have used E/fi = 2(1 + a) and = 2/(1 - a). The various contributions to the the integral 
are easily understood in terms of the different waves radiated into the cavity. Remember that ^ is the y 
wave number of the radiated waves in units of Ul- There are contributions to the integral for < < r 
corresponding to the radiation of transverse and longitudinal waves over all angles. In addition there is a 
contribution from the residue of the pole at -Fo(0 = which corresponds to the radiation of edge waves. 
Miraculously, the quantity in the square brackets numerically evaluates to 1.0 independent of ct in the allowed 
range —1 < a < 1/2, so that 

ro(w)=4fcu;. (151) 



3.3.2 In-plane Bending 

For an incident bending wave with unit displacement amplitude v in the y-direction there are two sources of 
radiation into the cavity: the oscillating moment 2^/2EIk^e^'^^* and the oscillating shear (tangential) force 
over the source width w in the cavity wall. The moment can be described in terms of the 
normal stress 2^/2Ek^e~^^*^y since Jy^dy = I. The tangential force has an additional power of A: '--^ ^/uj 
compared to this, but the radiation efficiency of the normal force is reduced by a power of kcW with kc the 
wave number of a propagating mode in the cavity kc ~ i-^/ct, since the two halves of the radiation source 
cancel at leading order. This means that the contribution of the normal stress to the power radiated is 
higher order in lu for small cj, and may be neglected. 



The analysis proceeds as in the previous section. We again use Eqs. (145) but now with = 
2V2E{I/w)k^e-''^^ and S„ ~ 0. This gives the equations for the mode amplitudes 

r^{ikTaT + ikLttL) + [r^ - 2){~ikTaT + iC^aL/kL) = (152a) 
*C(aT + ah) + {—ik'^aT/C + iC«i) = '^tw/p, (152b) 

which can be solved to yield the response to leading order 



(This result is analogous to Eq. (124) of Miller and Pursey ||l6|). This gives the average power radiated to 
leading order in small lo 

Prad ^ Re{-vY.;) = \j{2^/2EIk^/wf Re{iijjv /Y^t) (154) 
For the incident wave of unit amplitude we have 

{u, u', M, F) = (1, ik, Elk^, i£;/A:3)e'('=^-"*) (155) 
so that the average incident power is 

Pine = ^ Re [-iuj{Fu* + Mu'*)] = ojEIk^. (156) 



4(1 + .)^ r-^,^- 



The ratio gives the transmission coefficient 

%{oj)^\{wkf . 

3 [7r(l-cr) 7o FoKi) 

Again the quantity in the braces turns out to be unity, so we have 

1 



(157) 



%{uj) = ^{wkf. (158) 
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3.3.3 Flexural Modes 



The flexural displacement w{x,y) for a wave at frequency w satisfies the equation 



with = pd/Duj, cf. section 3.1.2. Expanding the cavity solution in transverse Fourier modes 



the Fourier amplitude w satisfies 



The solutions are w ^ e with 



Zn 



(159) 



(160) 



(161) 



(162) 



Solutions corresponding to a wave propagating away from the source at a; = or exponentially decaying to 
+CX3 are given by k± 



and 



k+ = + C2 



lyJC^ ~ C > K 

Thus the solution in the cavity can be written 

1 



The sources at a: = are a moment M{y) and an effective force per unit length V{y) 
M{y) = -D I + cr 



w-\F,{w/2)+F,{-w/2)) 



dx"^ dy^ 

ny) = -^UzT + (2"^)- 



dx^ 



dxdy^ 



(163) 
(164) 

(165) 

(166a) 
(166b) 



with Fourier transforms M(C) and V^(C)- The last two terms in the second equation are the corner forces. 
Matching these boundary conditions gives 



+ c\\ - a)\ w+iO - [K' - - a)] w^iO - -M{0/D, 
ik+ [K^ - C'(l - 't)] w+{0 - ik^ [K^ + C'(l - a)] w^O = -V{C)/D. 



The average power radiated is 



P 



w{x = 0, y)V{y) + 0{x = 0, y)M{y) dy 



(167a) 
(167b) 



(168) 



with 9 = —dw/dx the tilt angle and the dot denoting a time derivative. For oscillations e the average 
over time gives 



w{x = Q,y)V*{y)dy+ / e{x = Q,y)M* ijj) dy 



(169) 
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Evaluating the first integral in terms of Fourier expansions we find the power radiated by the force 



1 

Pw = -—Lulm 



(170) 



As in the scalar case, the imaginary part of this integral corresponds to the excitation of propagating waves 
for which ( < K and we may evaluate ^(C) to lowest nonzero order in Kw 



viC) 



w/2 



w/2 



dyV{y)e-'^y ~ - <V^i + • • • 



with 



^0= / dyViv), 



dyyV{y). 



(171) 

(172) 
(173) 



We only keep the second term for antisymmetric sources for which Vq is zero. Note that Vq is the total force 
normal to the plate, and Vi is the torque about the x axis, and these can be evaluated from macroscopic 
arguments. 



Similar arguments for Pg give 



Pe = — Lo Im 
47r 



with 



M(C) ~ Mo - <Mi + • • • 

with Mq, Ml the zeroth and first moments of M over the bridge end. 
Now we can calculate explicit results for the bending and torsion modes. 



(174) 



(175) 



Bending mode From Eq. ( |130| ) we see that an incident wave of unit displacement amplitude e^^''^ in 
the bridge gives oscillating sources on the edge of the cavity 



Mo = 2V2D(1 - a^)wk-e-' "e 



1^0 = 2\/2D(l - <7'^)wk"e-' 'e 



2 i7r/4 —iujt 

3 i7T/4 —iujt 



(176) 
(177) 



Note that we are using the macroscopic formulation to derive these expressions. It is somewhat subtle to 
directly use the expressions Eq. (166) since the y dependence of the mode structure cannot be ignored. We 
have verified that these expressions are reproduced using the long wavelength limit of the mode structure 
given by solving Eqs. (101,102). Now defining — (^/K and 

w±(C) = -2\/2we"/4w±(0, (178) 

(179a) 
(179b) 

(180) 



using the dispersion relation fc^ = iiT^/Vl — cr^, and then matching to the sources gives 

[1 +e{i- -)] ^+ - [1 - eii - <r)] = (1 - a^y/', 
-Vi+e [1 - e{i <j)] u+ ^Vl^ [1 + <^)] u^ = {i- o^f^ 

The power radiated is 

2 



P = uiDw'-k'-K\\ - a"-)- Im 

TT 



(m+ +U_)(1 - cr^) 



2^-1/4 _ / 
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Normalizing by the incident power Pi = ujDwk^il — ct^) gives the transmission coefficient 

Tq{uj) = kwli{a), 

where Ii is the integral 



h{a) = v/l-fi^-Im / {u+ + u-){l-a'')-^'^ + {^JTTen+-i^JT^u^) d(,. (182) 



For cr ~ 0.33, evaluating u± from Eq. (179) we find 

/i=2.3 ((7 = 0.33). 



(181) 



(183) 



Torsion Mode A unit amplitude mode 6 — e^^*^^ '^'^ gives the oscillating torque source 

T = 4Dw(l - cr)i/ce"''^*. (184) 



Here the wave number in the beam is given by (cf. Eq. (IIC)) 



kw — 



1 



2V6(1-'t) 



{Kw) 



(185) 



with ~ pdj Dlo as before. The torque Eq. (184) corresponds to a force source term given by a nonzero 
first moment Viof the tangential force 



7i ■ 



(186) 



This gives a source term on the right hand side of Eq. ( 167b| ) V{C) — ~^C^i ^^'^ the source term M in 
Eq. (167a) is zero. Now defining 



, 2u;Vl - cr 
w±{0 = ^= M±(?) 



(187) 



with = Q/K from Eq. (167) we find u± satisfies 

[1 + e^(l - a)] u+-[l- e^(l - cr)] ii_ = 0, (188a) 
- VTT^ [1 - ^2(1 - a)] u+ - i^l^ [1 + ^2(1 - a)] = e. (188b) 



Solving these equations for u± and so w±, substituting into Eq. (170) for the power radiated [Pg does not con- 
tribute for this mode) and normalizing by the incident power loDuP'K^^I — o'/2\/6 yields the transmission 
coefficient 



To(w) = kwl2{a) 



where I2 is the integral 



4 f°° 

/2(a) = (l-cr)-Im / {u+ + u^)^d^. 

71" J-00 



For a = 0.33, solving Eq. ( |188D for u± yields 

/2=0.6 (cr = 0.33). 



(189) 



(190) 



(191) 
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Mode 


ujik 0)/cTk 


uji/{ttct/w) 


To{k ^ 0) 




0), z = luw/ct 


Compression 


^2(1 + a) 


j 


■ikw 




Torsion 


2d/w 


8.014 fd\ 


l2kw 




z 


In- plane bend 


V ^ 


1 


^ikwf 




^^%3/2 


Flex-bend 


^kd 


2.886 rd\ 


Iikw 







Table 1: Long wavelength properties of the modes of a thin beam. See text for the details. 



4 Applications and Discussions 



The results for the long wavelength properties a thin plate beam of width w and thickness d are brought 
together in table |^. The second column gives the small k dispersion relation in terms of the propagation speed 
CT of the in-plane shear wave in a thin plate which is the same as the shear wave speed in the bulk medium. 
The third column gives the frequency cutoff ui of the lowest waveguide mode with the same transverse 
parity symmetry as the acoustic mode (this would be 2 A in the scalar model) . The fourth column expresses 
the small w, k energy transmission coefficient To of the acoustic mode in terms of the wave number k in 
the bridge and the width w of the bridge. This is useful in considering the Q of the fundamental vibration 
modes of the beam for which k is of order tt/L with L the length of the beam. The quantities Ii and I2 are 
Poisson ratio dependent numbers defined by Eqs. (182 , 190] ) and take on the values Ii = 2.3 and I2 = 0.6 
for a — 1/3. Finally, the fifth column re-expresses the small w, fc dependence of the energy transmission 
coefficient %, in terms of the frequency co. This form is particularly useful to estimate the reduction from 
the universal thermal conductance at low temperatures due to the strong scattering of the long wavelength 
modes by an abrupt junction. 



4.1 Implications for Heat Transport 

The low frequency asymptotic limit of the transmission coefhcient 7o{uj — s- 0) allows us to calculate the low 
temperature variation of the thermal conductivity. If we write the dependence as 

To{uj ^ 0) = A{ujw/ct)p (192) 

then relative to the universal low temperature one mode value K/T\^^ — 7r^k^/3h we have for each acoustic 
mode 

K/T , fkeTwY 3 n x^+^e" , 

where the prefactor A for each mode can be found from table The integral is just some p-dependent 
constant. Using the expressions in table [l] the result can be written in terms of the frequency uji cutoff of 
the first wave guide mode of the corresponding symmetry 



K/Tl V ^1 

The prefactor S is a numerical constant that depends on the Poisson ratio cr, but not on geometrical factors 
such as d/w. Since the thermal excitation of the waveguide modes occurs for fc^T > 0.2 x huji (cf. Fig. ^) 
this expression indicates to what degree the plateau in K/T becomes apparent as the temperature is lowered 
and the waveguide mode freezes out, before the reduced transmission coefficient at small frequencies begins 
to lower K/T to zero. The ideal low temperature universal value of K/T will be more evident for smaller 
powers p (cf. the comparison of the two scalar results for p ~ \ and p = 3 in Fig. (||)). This suggests that the 
compression and torsion modes will give contributions to K/T curves similar to the result for the stress free 
scalar model in Fig. H, without a well defined plateau at low temperatures (all have p = 1), and the in-plane 
bend mode (p = 3/2) will probably have no indication of a plateau. On the other hand for the flexural bend 
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mode {p ^ 1/2) the transmission coefficient increases more rapidly with increasing frequency %, ~ y^uj/uji. 
This will lead to a more rapid increase in K/T towards the universal value before significant excitation of 
additional modes occurs at T ^ huji/kB, leading to a more pronounced plateau. It should be noted that for 
this mode the plateau in K/T only develops at very low temperatures in the thin plate limit, reduced from 
the simple estimate hcT/kBW by the ratio d/w of the thickness to the width of the beam. 

4.2 Implications for Q 

Using the third column of the table and Eq. ( [lO| ) for Q we get the simple estimates Q L/w for the 
fundamental compression mode, torsion mode, and flexural-bending mode, and Q ^ (L/w)^ for the in-plane 
bending mode. Only for the in-plane bending mode is the isolation of the bridge modes from the supports 
sufficiently strong to give a large Q for accessible geometries (e.g. L/w < 100). The strong isolation of this 
mode is easy to understand physically: the wide supports are very rigid against bending motion in the plane. 

In many experiments on mesoscopic oscillators, modes other than the in-plane bending modes are used, and 
values of Q significantly higher than the value suggested by the geometric ration L/w are obtained. One 
way this is done is to use more complicated geometries, such as compound torsional oscillators arranged 
so that the amplitude of vibration in the bridge supports is reduced. Also, in oscillators at larger scales it 
is relatively easy to produce more rigid supports, for example by fabricating a bridge or cantilever making 
an abrupt junction to a three dimensional support, which can also be of a different elastic material, both 
of which will reduce the coupling to the support modes. In mesoscopic oscillators, where the geometry is 
typically etched out of a single material, and undercutting by the etch comprises attempts to make an abrupt 
junction to a three dimensional support, our estimates of the coupling will be more appropriate. 

Our estimates of Q suppose that all the energy communicated to the support modes is lost from the energy 
of the oscillator. This is not necessarily the case, for example if the support material is also of sufficiently low 
loss and isolated from the rest of the experiment. However, our results do show that when the bridge-support 
coupling is large, it is important to consider the dissipation properties of the support structures as well as 
the bridge, cantilever, or other oscillator that is the obvious focus of attention. 
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